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Graphene on mica probed by scanning tunneling microscopy: networks of ABA and 
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Graphene flakes are prepared on freshly cleaved mica by exfoliation in different environments and 
studied by scanning tunneling microscopy in ultra high vacuum down to the atomic length scale. On 
the monolayer, we observed an overall roughness of about 60 pm, after heating the sample to 180° C 
which is larger than on similarly prepared multilayers (24 pm) and larger than on monolayers which 
have been mapped recently by atomic force microscopy without heating (24 pm). A slight preference 
for a distinct wave length (A ^ 8 nm) of corrugations could be determined using an autocorrelation 
analysis. On few-layer graphene, a triangular network of partial dislocations separating ABC stacked 
and ABA stacked graphene was found similar to the networks occasionally visible on freshly cleaved 
HOPG. We found differences in the electronic structure of ABC and ABA stacked areas by scanning 
tunneling spectroscopy, i.e., a pronounced peak at 0.25 eV above the Fermi level exclusively in the 
ABA areas, which is largely responsible for the different apparent heights. 

PACS numbers: 61.72.Nn, 73.22.Pr, 61.48.Gh 



I. INTRODUCTION 

Since Geim and Novoselov succeeded in producing 
monolayer graphene flakes by exfoliation^, graphene de- 
veloped into a favourite research topic due to its two- 
dimensional nature and its unique combination of elec- 
tronic, mechanical and optical properties accompanied 
by chemical inertness.^ ^ 

Recently, trilayer graphene, which exhibits two differ- 
ent possibilities of stacking, came into focus^ ^, since 
it was predicted that an electric field E can open a 
band gap of up to 200 meV for ABC (rhombohedral) 
stacking, but not for the conventional ABA (Bernal) 
stacking.^Hi^ Interestingly, it has been shown, in addi- 
tion, that thicker layers of ABC stacking exhibit conduct- 
ing surface states in the absence of but an insulating 
bulk, very similar to weak topological insulators.l^In ad- 
dition, high-temperature superconductivity driven by the 
large density of states at the surface has been proposed 
for thicker ABC stacked layers. The two different stack- 
ings of trilayers could be distinguished experimentally 
by gate dependent infrared spectroscopy,^ by Raman 
spectroscopy,^^^^ and by differences in the temperature- 
or magnetic-field dependent transport properties.^ ^^^^^^ 
This revealed, e.g., that about 15 % of the sample 
area prepared by exfoliation is ABC stacked,'^ and that 
the predicted distinction in band-gap opening exi sts, al- 
t hough so far only with a band gap below 10 meVW^^^^ 
These methods, however, are limited in lateral resolution 
to about 1 /im, thus, they are not able to resolve finer 
stackin g o rders in graphene. Based on earlier work using 

scanning; tunneline; 

spectroscopy (STS) can distinguish ABA and ABC ar- 
eas exploiting a peak in scanning tunneling spectroscopy 
found only on the ABA stacking areas. We identified a 



region on multilayer graphene on muscovite mica which 
exhibits a triangular network of dislocation lines separat- 
ing ABC and ABA areas with sizes down to (200 nm)^. 
For such a network, one can anticipate, based on the the- 
oretical predictions, that a vertical electric field would 
result in a triangular network of graphene quantum dots 
with tunable barriers. 

Secondly, we investigated the corrugation of graphene 
on mica prepared with a heating step to 180-190° C. 
The rms value of the surface roughness was about 60 
pm for monolayer and 25 pm for multilayer and a slight 
preference for a wavelength of 8 nm in the corruga- 
tion pattern was found. The roughness of the mono- 
layer was larger than the 24 pm determined previously 
by atomic force microscopy (AFM) without heating the 
sample prior to the measurement. Water intercalation 
between graphene and mica, leading, in addition, to a 
p-type doping of graphene, has been regularly found 
by AFM after preparation in air or even at strongly re- 
duced humidity. ^^ ^^ We tried to reduce the density of 
water below the graphene by exfoliating the graphene in 
Ar atmosphere with less than 2ppm of water, however, 
without success. The water below the graphene remains 
to cover about 50 % of the area below graphene which 
implies that the water is not trapped during the exfolia- 
tion. 



II. EXPERIMENT 

Graphene samples were prepared by mechanical 
exfoliation^ on freshly cleaved muscovite mica (high- 
est quality "VI", Piano GmbH, Wetzlar, Ger- 
many), an aluminosilicate with the formula unit 
KAl2(Si3,Al)Oio(OH)2^^ and a layered structure consist- 
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ing of one layer of octahedrally coordinated Al^+ions sur- 
rounded by two layers of tetrahedral Si^+ ions. Every 
fourth Si^+ ion is replaced by an Al^+, resulting in an 
excess negative charge, which is compensated by K+ ions 
connecting the triple layers. Mica can be cleaved along 
the K+ layers with large areas (several 100 /im^) of mono- 
atomic flatness. However, the density distribution of K+ 
on the two cleaved surfaces is not perfectly homogeneous 
leading to hydrophilic properties which trigger the for- 
mation of potassium car bonat e islands on the surface, if 
ambient water is present .^^^^ 



In order to reduce the influence of water, we cleaved 
mica and deposited graphene on top by exfoliation not 
only under ambient conditions, but also within a dry air 
box with reduced though not explicitly controlled water 
density as well as within an Ar box with a measured H2O 
concentration below 2 ppm. 



Thin graphene flakes were identified afterwards using 
optical polarization microscopy under ambient conditions 
as described by Dorn et. al..^^ The number of layers was 
determined by Raman spectroscopy at laser excitation 
wavelength of 532 nm.^^ To avoid extensive heating of 
the sample, the laser power was kept below 2 mW and the 
measurement time was 5-10 s. To determine the number 
of layers in the case of few-layer flakes as presented in 
section |IIIB[ AFM was used to measure the height dif- 
ference between a monolayer area of graphene and the 
multilayer area of interest. 



To avoid possible chemical residues on the graphene 
surface resulting from lithography,^^ indium micro- 
soldering was used for contacts as initially described by 
Girit et. al-.f^This required heating of the sample to 180- 
190° C in order to melt the In onto the graphene surface 
properly, which might have influe nces on the observed 
corrugation as discussed in section III Cp ^ 



The STM measurements were carried out at 300 K 
in a home-built UHV system operating at a background 
pressure of 1-2-10"^ Pa. Tungsten tips were prepared by 
electrochemical etching and subsequently short heating 
inside the ultrahigh vacuum chamber to remove the tung- 
sten trioxide layer formed during the etching process.SSl 
The positioning of these tips above the graphene flakes is 
achieved using a long-range optical microscope with reso- 
lution down to 10 /im. A bias voltage V ofV = 0.2- 0.6 V 
orV = — 0.2 — 0.5V was applied to the graphene sample 
and currents of / = 0.05 - 0.5 nA were used for constant- 
current images. The dl /dV-cuiYes were acquired by 
applying an additional ac-voltage with a frequency of 
1.4318 kHz and a modulation amplitude of Vmod = 30 mV 
and detecting the resulting in-phase ac current by a lock- 
in amplifier. 
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FIG. 1. (a) + (c) Optical images of samples prepared by exfo- 
liation in the dry box recorded with polarization microscopy; 
layer thickness n in monolayers of graphene is marked as nL. 
(b) + (d) Corresponding Raman spectra of the regions marked 
in the optical images as IL (Monolayer), 2L (Bilayer) and 5L 
(Multilayer). 



III. RESULTS AND DISCUSSION 

A. Raman spectroscopy and atomic force 
microscopy 

Figure [l] shows optical images of sorae of the prepared 
samples along with the corresponding Raman spectra. 
The number of layers is identified as usual by analyz- 
ing the shape and position of the 2D peak as described, 
e.g., by Malard et. al..^^ For thicker layers the steps be- 
tween monolayer and adjacent layers are imaged by tap- 
ping mode AFM in order to determine the height as given 
in Fig. [1] Fig. |2] shows a tapping-mode AFM image of 
the sample prepared in the argon box at relative humid- 
ity below 2 ppm. The monolayer is folded back at the 
edge and several first-layer and second-layer water ad- 
layers as well as larger bubbles are visible, which are be- 
lieved to be trapped between the graphene layer and the 
mica substrate. 3i,4i | rpj^^ ^^^^ covered by water under- 
neath the flake is roughly 53%, i.e. even larger than the 
area obtained after preparation under ambient conditions 
(22%). Like Shim et. al.^^, we were unable to see any re- 
lationship between relative humidity and water coverage 
during exfoliation. As w ater coverage can be even larger 
than a single layer,'^^^ the AFM images do not exclude 
the possibility of more than one water adlayer below the 
graphene directly. However, the distinct topography be- 
tween the lower and the upper areas of graphene disfavor 
such an assumption. A similar distinct topography, in 
particular, a distinct apparent height of adjacent steps, 
is visible in STM measurements (s. Fig. [3|. Thus, most 
likely, a submonolayer of water is trapped between the 
graphene and the mica even after preparation in nearly 
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FIG. 4. Sketch of ABC stacked trilayer graphene with the 
three layers colored differently as marked. Only the carbon- 
carbon bonds are shown as lines. The arrow labeled 'a' 
marks a Burgers vector of a possible dislocation in the up- 
per graphene layer separating two areas of the same stacking. 
This dislocation can be separated into two shorter (and thus 
energetically more favorable) partial dislocations (labeled 'si' 
and 'S2'). Since 'si', for example, transforms an A-layer into 
a C-layer, adjacent stacking- types of ABC and CBC result, 
i.e. an area of different stacking. 



FIG. 2. Tapping mode AFM image acquired under ambi- 
ent conditions of a graphene monolayer on mica with folded 
edge on the left; cantilever stiffness: ^ 40^, frequency: 
300.25 kHz, set-point: 52%, free amplitude: 31 nm; the sam- 
ple was prepared in an argon box; water underneath the 
monolayer (elevated areas) covers about 53% of the area be- 
low the graphene. 




FIG. 3. STM image showing water underneath a monolayer 
graphene (s. Fig.[T]^a)) prepared by exfoliation in the dry box. 
V -0.5 V / = 0.2 nA. 



H2O free environment. This is surprising since the water 
adlayer on uncovered mica has been shown to disappear 
below a relative humidity of 20 We do not know the 
origin of this water. Possibilities could be that the water 
diffuses below the graphene after preparation, when ex- 
posed to ambient conditions again or it diffuses out of the 



mica, e.g., due to the pressure applied during exfoliation. 



B. Stacking- faults in few- layer graphene on mica 

For more than two layers, the stacking order of 
graphene has to be taken into account. There are several 
ways of stacking. The most stable one is called Bernal- 
stacking showing an ABA sequence as illustrated in Fig. 
|4] Another obvious possibility is rhombohedral stacking 
with ABC sequence being energetically slightly less favor- 
able by about 0.18 meV per atom.^^ For both stackings, 
the distances between carbon atoms within the graphene 
sheet (1.418 A) and the separations of graphene layers 
(3.348 A) are identical.SSl 

The lateral energy barrier between these two stack- 
ings has been calculated to be about 2 meV pe r atom 
only, which favors a network of stacking faults.^ There- 
fore, partial dislocations with a Burgers vector si = 
|a<1100> are required to transform ABA graphite into 
ABC graphite and back again (s. Fig. [4|. A network 
of partial dislocations connecting triangular ABA and 
ABC areas has indeed been observed in highly oriented 
pyrolytic graphite (HOPG) by electron microscopy.'^^'fl^l 
The partial dislocation character has been proven di- 
rectly in dark-field images. STM images of such dis- 
location networks were recorded subsequently revealing 
an apparent height contrast between differently stacked 
areas.l^^K^ Kobayashi et. al.^^ showed, in addition, that 
this contrast changed in intensity across a buried step 
edge, which was qualitatively explained by a difference 
in the local density of states (LDOS) resulting from an 
interference of electronic waves normal to the surface. 
STM images of the few-layer graphene area shown opti- 
cally in Fig.[l](a) are displayed in Fig. [5] revealing an area 
with a triangular pattern. Triangles pointing in differ- 
ent directions exhibit apparent heights differing by 1-4 A. 
Due to the similarity of theses patterns to the patterns 
found previously on HOPG,^^ 44 47 assign it to a 



line defect 

/ 



4 




FIG. 5. STM images of dislocation network on few- layer (6 ± 1 layers) area of graphene on mica, (a), (b) are recorded on the 
same area with 2 hours in between; V — 0.4 V, / = 0.1 nA, scan speed: 2 white circles highlight a defect which marks the 
same area in (a) and (b); a line defect is indicated in both images; notice the change of the network towards a more regular 
structure from (a) to (b); (c) Zoom into the lower left edge of (b); white lines mark dislocation lines, which according to the 
indicated curvature presumably encircle ABA and ABC areas as marked; (d), (e) are recorded subsequently on the same area 
with different voltage polarity; white circles mark a defect, i.e. the same position in both images; (d) V = 0.3 V; / = 0.15 nA; 1 
i^; (e) V = -0.3 V; / = 0.15nA; 1 notice the contrast inversion between (d) and (e). (f) STM image recorded at V = 0.4 
V; / = 0.2 nA; 0.6 ^; during scanning, a 5 V pulse is applied at the white line as marked, leading to contrast inversion. 



dislocation network separating triangular areas of ABA 
and ABC stacked graphene. Fig. [s] (a) and (b) show 
two large-scale images, which have been measured subse- 
quently on the same area of the sample. Obvious changes 
in the pattern towards a more regular network emphasize 
the mobile nature of the dislocation network, which is a 
well-known property of partial dislocations.^^ Within the 
first image, some dislocation lines appear to be pinned to 
the defect line marked in Fig.[5](a), but are depinned two 
hours later in Fig. [s] (b). We checked carefully that the 
appearance of such a relaxed structure does not depend 
on scan direction, which, however, does not exclude that 
the depinning was initiated by the scanning process. The 
large scale mobility of the network underlines that the 
apparent height is not caused by a topographic feature, 
which is unlikely to be moved strongly at 300 K once and 
then remain stable. The triangular shape of the pattern 
is, in addition, in line with the three possible directions 
of Burgers vectors transforming ABA into ABC areas. 
In Fig. [5] (c), the shape of a dark and a bright region 
is marked by white lines. One observes that the bright 



one is a triangle with convex edges, while the dark one 
exhibits concave edges. Based on this shape, the stack- 
ing order of the different areas can be assigned. Since 
ABC graphite is energetically less favorable, there is 
a force on the dislocation lines in the direction of ABC 
stacking, resulting in concave shapes for these areas. In 
turn, this identifies the darker areas observed at positive 
sample bias as ABC regions or, more generally, as the 
stacking fault regions. The crossing points of the dislo- 
cation lines, called nodes, require an energetically costly 
local AA stacking or a vacancy due to symmetry. These 
nodes often appear brighter than the surrounding area 
(e.g Fig.|5](d)). 

The contrast between ABA and ABC stacked areas as 
shown in Fig.jsjd) and (e), can invert with bias polarity. 
This was not observed on HOPG, where only a gradual 
change was found. Moreover, the contrast could even 
be inverted by voltage pulses changing the shape of the 
tip. This is shown in Fig.jsjf) and verified by imaging the 
complete area prior and after the voltage pulse with the 
same tunneling voltage of 0.4 V (not shown). Since the 
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FIG. 6. dl /dV-specii a, acquired on ABA-stacked (solid) and 
ABC-stacked (dashed) regions. Spectra of the same color were 
acquired on adjacent regions directly after each other. The 
black, red and blue measurements are stabilized at 1/ = 1 V 
and / = 0.2 nA, the pink ones at V = 0.5 V and / = 0.2 nA, 
Vmod = 30 meV. Measurements are shifted vertically for clar- 
ity. Vertical lines mark minima. Notice that the pink spectra 
are slightly shifted upwards with respect to the other ones. 
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FIG. 7. (a) Small scale STM image of multilayer graphene 
with stacking faults, V = 0.25 V, / = 0.15 nA; (b) same area 
as in (a), V = 0.6 V; / = 0.15 nA. (c) Voltage dependency of 
the averaged height difference between ABA and ABC areas 
Az and of the average height of the nodes with respect to the 
highest stacking area in the surrounding (nodes); / = 0.15 nA. 
(d) Current dependency of the same quantities as displayed 
in (c); V = 0.4V. 



tip where we occasionally observed the contrast change 
by pulsing was rather unstable as visible in the STM im- 
age and exhibited an unusually large decay length in 7(2:) 
curves, we believe that the tip chemistry and, thus, the 
surface potential of the tip changed significantly during 
the pulse. Importantly, all these observations strongly 
suggest that the observed contrast is not of topographic 
origin but of an electronic one. 

To further analyze the electronic differences, dl/dV- 
spectra, which are proportional to the local density of 
states (LDOS), were acquired on ABC and ABA areas 
as shown in Fig. [6j The spectra exhibit several min- 
ima marked by vertical lines and an obvious peak at 
V = 0.25 V only on the ABA areas. Also the positions 
of the minima are different for differently stacked areas 
as marked. Since the curves are acquired alternately on 
ABA- and ABC-stacked areas, one can exclude that these 
spectroscopic differences are caused by a tip change. The 
global minimum is found around mV for ABA areas and 
at 125 mV for ABC areas. This might indicate a relative 
offset of the Dirac point energy £^d- The peaks and dips 
are slightly moved upwards by about 20 mV in the spec- 
tra acquired at smaller tip-surface distance (pink curves) . 
This indicates an upwards band bending at the graphene 
surface induced by the work function mismatch between 
tip and graphene and implies an electric field penetration 
into the graphene sample. No band gap is observed nei- 
ther on the ABA nor on the ABC areas, which means that 
a band gap, if present, is smaller than about 100 meV be- 
ing the energy resolution of the experiment. Theoretical 
calculations^ ^'^^ ' ^^* of ABA and ABC stacked graphene 
find complex band structures changing with electric field. 



which, however, are largely symmetric with respect to 
£^D- We checked that they cannot explain a single peak 
at +0.25 eV, which is only found in ABA areas, using 
reasonable assumptions for the electric field penetration 
depending on tip-sample distance and screening. Since 
the peak is observed close to the tentative £^d of the 
ABC area, i.e. where a gap exists in ABC stacked ar- 
eas, it might indicate a confinement of electronic states 
within the ABA area (see below). 

Based on the spectroscopic data, we can show more 
quantitatively that the apparent height difference be- 
tween ABA and ABC stacked areas is of electronic ori- 
gin. Therefore, we first measure the dependency of Az 
on sample voltage V and current /. Figure [Tj^a) and (b) 
show part of the dislocation network imaged at two dif- 
ferent V ^ where the height difference between the darker 
areas (ABC) on the left and the brighter areas (ABA) on 
the right has been used to measure Az. The results are 
displayed in Fig. [tJ^c) and (d). In addition, the height 
of the nodal points, is determined with respect to the 
surrounding area appearing the highest and is shown as 
'nodes' in Fig. [tJ^c) and (d), too. Obviously, the current 
dependency of Az is weak while the voltage dependency 
is pronounced. The contrast between ABA and ABC ar- 
eas inverts with voltage polarity, increases with decreas- 
ing |V|, and is largest around V — 0.25 V. Similar trends 
are observed for the contrast of the nodes, except that 
the nodes always appear higher than the surrounding. 
Partly, the dislocation lines exhibit a distinct contrast as 
visible in Fig. [7|^b), but we did not evaluate this contrast 
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in more detail, because it appears only at certain V. 

Next, we calculate the apparent height difference ex- 
pected from the dl/dV curves. According to Tersoff and 
Hamann^^, the tunneling current / at voltage V depends 
on the LDOS of the tip p\E) and the LDOS of the sam- 
ple at the position below the tip p^{E) according to: 

/oc / p'{EF^e)-p'{EF-eV^e)-T{e,V,z)de, (1) 



Here, Ep is the Fermi level, z is the distance between the 
tip and the sample surface, and T(e, V, z) is the transmis- 
sion coefficient of the tunneling process being exponen- 
tially dependent on z. Since a constant / is used to sta- 
bilize the vertical tip position with respect to the sample, 
a difference in LDOS of the different sample areas results 
in a change of z and consequently appears as a height dif- 
ference Az in constant-current images, even if the sample 
is flat. To estimate the expected Az at voltage the 
dl/dV curves are integrated between V and the applied 
voltage Vmeas ^ud translated into a height difference us- 
ing /(z)-spectra. The procedure is illustrated in Fig. [s] 
exemplary for the pink curves of Fig. [6j As can be seen 
by comparing Fig.jsj^a) and (b), the dl/dV values in the 
completely filled area of the curve are larger for ABA- 
stacked graphene than for ABC-stacked graphene. Thus, 
the resulting height difference must increase by reducing 
the voltage from V^tab to Vmeas- This is indeed observed 
experimentally. To determine the expected Az quan- 
titatively, the filled area underneath the di/dV- curves, 
ranging from V = OV to Vmeas ^ is divided firstly by the 
hatched area, ranging from V = OV to Fgtab- Then, the 
resulting ratios obtained for ABA and ABC area are di- 
vided by each other and are used to determine a height 
difference Az from the exponential I{z) curve as illus- 
trated in Fig. [sj^c). For the V^meas in the example, an 
additional height difference of 111pm results, which has 
to be added to the height difference observed at Vst^h to 
determine Az. Notice that the exponential I{z) curve 
guarantees an independence of the found Az from the 
chosen z value in the /(z) curve. 

Figure [9] shows the result of such a calculation as 
a function of Vmeas using dl/dV curves stabilized at 
V'stab = 1 V (red line) in comparison with the determined 
Az from Fig. 7|^c) (black squares). Dashed lines mark 
the error bar determined from five pairs of dl/dV curves 
recorded on different positions. As the topographic data 
at y = 1 V became rather noisy, the height difference 
dit V = 0.6 V was added as the offset to determine the 
calculated heights. Note that the simulated height differ- 
ence changes very little between V = 0.6 and V = IV. 
The characteristics of the measured Az{V) curve, i.e., 
the change of contrast around V = V and the decrease 
of contrast with increasing |V|, are well reproduced by 
the calculation, although the experimental |A2:| values 
are slightly larger. The reasonable agreement excludes 
a dominating influence of a mechanical interaction be- 
tween the tip and the graphene as well as of the elec- 
tric field of the tip onto the apparent height, since both 



are changing with distance being significantly different in 
the constant-current image used to determine the experi- 
mental Az and in the dl/dV curve used to determine the 
calculated Az. Since the contrast inversion with voltage 
polarity excludes a dominating influence of topography, 
we conclude that the different LDOS is primarily respon- 
sible for the contrast between ABA and ABC areas in 
constant-current images. 

The question remains: what causes the differences in 
LDOS, respectively the LDOS peaks? Intrinsic differ- 
ences of the band structure could be excluded as dis- 
cussed above. A Fabry-Perot type interference, similar 
to a proposal by Kobayashi et. al.^^, was considered. 
The idea is the following: electrons tunneling into the 
graphene multi-layer are reflected back at the interface 
between graphene and mica forming an interference pat- 
tern with the incoming electrons. Depending on the en- 
ergy of the electrons, this interference can be constructive 
or destructive causing peaks and minima in the LDOS. 
The energy of these peaks is straightforwardly calculated 
using the energy dispersion of graphite orthogonal to the 
surface: 
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with m±_ = 2c^j^ 7 the interlayer distance c = 335 pm and 
the interlayer hopping parameter 71 = 0.39 eV as well as 
the condition for constructive interference: 
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with n G N and N being the number of layers in the 
graphene sheet. 

The stacking fault might scatter the electrons isotrop- 
ically due to the band gap in its interior prohibiting the 
vertical interference, and, thus, explaining the absence 
of strong peaks in the ABC region. However, for the 
determined layer thickness = 6, the four lowest en- 
ergy peaks would be at about -0.7 eV, -0.35 eV, +0.2 eV, 
+0.9 eV with respect to in clear disagreement with 
the measured spectra exhibiting a single strong peak at 
about 0.25 eV and maybe a smaller one at 0.4 eV. Layer 
numbers N between 5 and 7 have been checked, in accor- 
dance with the maximum uncertainty of the layer thick- 
ness determined by AFM, but without success. Thus, we 
do not have any quantitative model for the strong peak, 
albeit it could be that an additional lateral confinement 
by the gap within the interior of the ABC area strength- 
ens the vertical Fabry-Perot resonance in the ABA region 
which is closest to £^d, i-e. the one at 0.2 eV. 



C. Corrugations of graphene on mica 

Finally, we consider the corrugation of graphene on 
mica. Even freely suspended graphene was found to be 
corrugated by about 100 pm using electron microscopy.!^ 
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FIG. 8. Illustration how to determine the expected height difference between ABA- and ABC-stacked areas from dl /dV-spectia; 
(a),(b) dl/dV -spectra: of ABA- and ABC-stacked areas as marked (same as pink curves in Fig.|6|; the hatched area up to Vstah 
has to be integrated to determine the current relevant for the stabilization during the dl/dV curve, while the completely filled 
area up to Vmeas determines the current at Vciieas; the ratio between filled and hatched areas is given above the curves; (c) 
sketch of the distance dependence of the current; the ratio of the ratios from (a) and (b), 0.7, is marked symbolically and at 
the /-axis; it is translated to a height difference by the I{z) curve as marked. 
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FIG. 9. Black Squares: height difference /\z between ABA- 
and ABC-stacked areas measured in constant- current images 
(same data as in Fig.[7|^c)); black solid lines serve as guide-to- 
the-eye; red solid line: averaged height differences calculated 
from five dl / dV-speciid. stabilized at Vst&h = IV according 
to the method illustrated in Fig. [Sj red dashed lines: error 
bars resulting from fluctuations for dl / dV-speciid. taken at 
different positions on the sample. 



This has firstly been explained by a coupling between in- 
plane and out-of-plane phonons,-^^ stabilizing graphene 
against thermal vibrations, that would normally crumble 
a two-dimensional material.^ Later experiments showed 
a strong influence of the preparation conditions on rip- 
pling, in particular, of the heating process. To an- 
alyze the influence of substrates, the graphene mor- 
phology has be en m easured, e.g., for gra phene de- 
posited on Si02 P^ * ^^^ ^ muscovite mica,^!^ hexagonal 
boron nitride,'^ and graphene grown on copper foilsP^^^^ 
Graphene on copper exhibited very large roughness due 
to the roughness of the underlying copper substrate and 
additional strain which builds up during cooling after 



growth and gets released by the formation of ripples and 
wrinkles. Exfoliated graphene on silicon dioxide also 
exhibits a large roughness of drms = 140 — 320 pm due to 
the underlying substrate and partly an intrinsic rippling 
observed, if thegraphene is not in complete contact with 
the 8102 ^^^^^ ^ ^ In the latter case, the suspended ar- 
eas could be hysteretically lifted by electrostatic forces.^ 
Graphene on mica was studied firstly by AFM revealing 
a roughness of 24 pm for the monolayer, which, however, 
might be limited by the lateral resolution of the AFM 
given as 7 nm.l^ STM data of graphene on mica have also 
been published with an intentionally complete interca- 
lated water layer between the mica and graphene making 
a roughness analysis difficult. Nevertheless, a preferential 
roughness value of 15 pm has been gwenW^ The small- 
est graphene roughness was observed so far on the flat 
van-der-Waals surface of hexagonal boron nitride (0001) 
being between 2 and 17 pm, whic h inte restingly, led to 
reduced potential fluctuations, too.l^^^^ 

Fig. [To] (a) and (b) show typical STM images of mono- 
layer and few-layer graphene on mica recorded in ar- 
eas neither exhibiting signs of underlying water nor of 
a stacking fault. A line fit and a low-pass filter in the 
Fourier transformation with a cut-off wavelength of 3nm 
and a first-order butterworth profile were applied to re- 
move the influence of occasional changes in the tip and 
high-frequency noise, respectively. Atomically resolved 
images, also with a low-pass filter applied, are displayed 
as insets exhibiting the hexagonal and triangular symme- 
try as expected for monolayer and multilayer graphene. 
Triangular symmetries, however, are partly observed on 
the monolayer, while hexagonal rings are seen on the 
multilayer, too, as shown in Fig. 11 This renders an 
assignment just by the appearance of the atomic lattice 
ambiguous j ^^ * ^^ * and calls for the additional characteri- 
zation by Raman spectroscopy as shown in Fig. [l] The 
height profiles shown in Fig. 10 ^c) reveal that the cor- 
rugation on the monolayer is significantly larger than on 
the multilayer. The rms roughness cr can be deduced 
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FIG. 10. (a) STM image of monolayer graphene on mica; y = 0.4 V, / = 0.2 nA, cr = 53.6 pm; (b) STM image of few-layer 
graphene on mica; V = 0.3 V, / = 0.15 nA, cr = 22.2 pm; insets in (a) and (b): atomically resolved images; (a) V = —0.4 V, 
/ = 0.3 nA, (b) V = 0.4 V, / = 0.1 nA; (c) height profiles along the lines marked in (a) and (b) in corresponding color; (d),(e) 
autocorrelation images AC{x,y) of (a), (b), respectively without any filtering beyond a line fit; (f) radial averages of (d) (red 
line) and (e) (blue line) with preferential wavelength marked; (g) histograms of the height values found in (a) (red line) and 
(b) (blue line) normalized to the average height. 




FIG. 11. Atomically resolved STM images on the monolayer 
(a) and the multilayer (b); different symmetries are indicated 
as black hexagons and triangles in both images. Images were 
recorded at (a) V = -0.4 V, I = 0.3 nA; (b) y = 0.4V,I = 
0.1 nA. 

from the corresponding histograms shown in Fig. [lOjg). 
For few-layer graphene, we find a = 24 ± 2 pm by av- 
eraging several images, which is identical to the value 
found on monolayer graphene by AFM without any heat- 
ing (24 pm)! ^^^ The monolayer graphene in our STM 
measurements, however, exhibits a larger roughness of 
a = (61 ± 13) pm, where the error bar indicates fluctua- 
tions between different measurements. The larger rough- 
ness could be due to the fact that our samples, unlike 



the samples of Lui et. al.,!^ were heated to 180-190°C 
during contacting which can lead to additional corru- 
gations induced by the different thermal expansion of 
graphene and substrate.^^ The thermal expansion coeffi- 
cient (TEC) of mica is 7.47 • 10" ^K,^^ while the TEC of 
monolayer graphene deduced from Raman spectroscopy 
is about -4-10-^/K between 0-200°C.^^ The latter value 
is negative and larger than the TEC of graphite, which is 
-1.4 • IQ-VK at 0°C and -0.8 • IQ-VK at 200°C.'^o The 
heating, thus, leads to a tensile strain of e = 0.2 % in the 
monolayer and e = 0.16 % in the multilayer of graphene 
at T = 180° C, if we assume a relaxed graphene layer af- 
ter transfer at 300 K. If the strain is relaxed at the high 
temperature, a compressive strain with the same absolute 
value of e builds up during cooling, which might induce 
additional corrugation. Using a triangular approxima- 
tion of the rippling with periodicity A and corrugation 
height /i, one gets: 

Deducing from Fig. [l0|c) that the dominating period- 
icity is A 8 nm , we find h o:^ 250 pm in good cor- 
respondence with the height of the corrugations on the 
monolayer visible, e.g., in Fig. [Tojc). Thus, most likely 
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TABLE I. Comparison of rms roughness a for mono- and few- 
layer graphene A on mica and Si02 



System 


roughness rms / pm 


monolayer / mica 


61 ± 13 (this work) 




2429 


few- layer / mica 


23 ± 2 (this work) 


monolayer/ Si02 


140_32d^^l^ 



the additional corrugation is induced by heating during 
the preparation process. It is of course possible that the 
lateral resolution of 7nm in the AFM measurements of 
Lui et. al.,^ slightly decreased the observed corrugation 
with respect to the real one. However, we simulated the 
topography of the image shown in Fig. [To] (a) as if imaged 
with a spherical tip always touching the surface but not 
deforming it and resulting in the apparent width of a step 
edge given in Ref.^^. However, the roughness remains the 
same (53 nm). This verifies that the roughness is indeed 
larger by a factor of three after our sample preparation. 

The observed change in corrugation height between 
monolayer and multilayer by a factor of three is simi- 
larly observed for graphe ne on Si02, albeit with larger 
a values as given in table This relation is probably 
caused by the larger stiffness of the multilayer. 
Finally, Fig. [lOjd) and (e) show the autocorrelation im- 
ages T/jof Fig. flOja) and (b) calculated as: 



AC{x,y) 



M-l- 
k=0 



''Y.t'o-' z{k,l)xz{k + x,l + y) 



{M-l-x)-(N-l-y) 



(5) 

where k and / label the x- and ^-positions of individual 
pixels up to the size of the STM image M, N and z{k, I) 
is the measured topographic height at the position (/c, 
Prior to determining AC{x^ y)^ we applied a line fit to the 
STM data in order to remove the influence of occasional 
changes of the tip, but we checked carefully that this line 
fit did not induce any novel anisotropy into the AC{x^ y) 
images. No further correction or filtering was applied to 
the raw STM data used to determine AC{x^y). 

Figure [lOj^d) and (e) show AC{x^y) for monolayer and 
multilayer, respectively. The AC{x^y) images are cut 
down to the central (200 x 200) nm^, since they become 
less meaningful at larger (x^y) as the number of pixels 
contributing to Eq. [5] is reduced. Both AC-images ex- 
hibit a complex pattern with certain regularities. The 
smallest regular scale can be determined by radial av- 
eraging the AC{x,y) data leading to the curves AC{r) 
with r = \J x^ ^ y^ as shown in Fig. 10 'f). The most ob- 



' x^ y^ as shown in Fig. 
vious oscillation occurs with a wavelength of 6 nm on the 
monolayer. This number slightly varies for different ar- 



eas (10 images evaluated) between 6nm and 10 nm. Such 
a preferential wavelength corresponds to the preferential 
distance between two hills or valleys of the corrugation. 
The amplitude of this oscillation is only 0.3 % of the value 
AC(0 nm,0 nm) = a^, which would imply that this pref- 
erential wavelength contributes only with a tiny fraction 
(5-6 %) to the total corrugation. A more detailed anal- 
ysis showed that phase differences of the oscillation in 
different directions considerably reduce the amplitude in 
AC(r), i.e. individual profile lines of Fig, p^d) revealed 
about the same oscillation period with a larger amplitude 
of 3-5 % implying that, at least, 15-20 % of the corru- 
gation happens at this wave length. Indeed, the red line 
section of the STM image of the monolayer in Fig. 10 'c) 
exhibits a preferential distance between peaks of 6-8 nm 
directly. Of course, we checked that the preferential wave 
length is not induced by a noise frequency of the instru- 
ment and that such a wavelength is not found after using 
our programme to determine AC(t^ for an image with 
random values on a grid of the same number of pixels. 
Thus, we have to conclude that preferential wavelengths 
develop during preparation similar as already observed 
for graphene on Si02!^ 



IV. CONCLUSION 

STM measurements on monolayer and multilayer 
graphene on muscovite mica are presented. On the mono- 
layer, we observe a short-scale corrugation with a prefer- 
ential wave length of 8 ± 2 nm and an rms roughness of 
63 ±13 pm, the latter being larger than determined previ- 
ously by AFM.2^ The roughness is reduced on few-layer 
graphene by about a factor of three. On a particular 
area of the multilayer, a regular pattern of triangles is 
attributed to a network of partial dislocations separat- 
ing areas of ABA and ABC stacking. STS measurements 
revealed that the STM contrast between these areas, ex- 
hibiting a size of about (200 nm)^, is of electronic origin 
largely dominated by a peak in the local density of states 
around 0.25 eV above the Dirac point. This peak is only 
found on ABA areas and might be related to a confine- 
ment effect resulting from the topological character of 
ABC graphene. 
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